Roughness exponent in the fracture of fibrous materials 
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PACS 



In this paper, a computational model in (2+l)-dimensions which simulates the rupture process 
of a fibrous material submitted to a constant force F, is analyzed. The roughness exponent £ at the 
boundary that separates two failure regimes, catastrophic and slowly shredding, is evaluated. In the 
catastrophic (dynamic) regime the initial strain creates a crack which percolates rapidly through 
the material. In the slowly shredding (quasi-static) regime several cracks of small size appear in all 
parts of the material, the rupture process is slow and any single crack percolates the sample. At 
the boundary between these two regimes, we obtained a value £ ~ 0.42 ± 0.02 for the roughness 
exponent, in agreement with results provided by other simulations in three dimension. Also, at this 
boundary we observed a power law behavior on the number of cracks versus its size. 
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I. INTRODUCTION 



The fracture process in disordered materials is a sub- 
ject of intensive research and has attracted much scien- 
tific and industrial interest The fracture process 
is extremely sensitive to disorder, which may act as ob- 
stacles difficulting the propagation of cracks through the 
material. Thus, disorder has a strong influence on the 
roughness of the fracture surface. Several computational 
models have been constructed to study the phenomenon 
of fracture of these materials, such as, the extensively 
studied fuse Q] and the well-known fiber bundle models 
, created from the pioneer work of Daniels |l0| . 

Experiments have shown that the fracture surface in 
disordered materials can often be described by self-affinc 
scaling In this case, the roughness of the frac- 

ture profile can be characterized by a roughness expo- 
nent £. Some experimental works have claimed that the 
roughness exponent £ has a universal value of 0.8 [[hi] |l6| . 
However this universality was questioned by Milman et 
al [0, which experimentally found a roughness expo- 
nent closer to 0.5. From the theoretical point of view, 
numerical models have been searched in order to evaluate 
the roughness exponent £. Simulations have shown that 
£ ~ 0.7 in two dimensions |l^,[l£| and that C ranges from 
0.4 to 0.5 in three dimensions [2^,^lJ. Nowadays there 
is a conjecture relating the smaller and the higher value 
of the roughness exponent £ to the speed of crack prop- 
agation through the sample pl]j . The greater value has 
been associated with a high speed of crack propagation 
and interpreted as a dynamic regime. In contrast, the 
smaller value of the exponent was related with a quasi- 
static regime, where the dynamic effects of the propaga- 
tion are negligible. Experimentally, the smaller and the 
higher value of £ are connected to the length scale at 
which the crack is examined. The smaller value is as- 
sociated to small length scales and the higher value is 
connected to large length scales Several models 

used to study the characterization of the fracture surface 
by the roughness exponent do not consider the influence 



of the temperature on the fracture process. In this paper 
we present a fibrous model in (2+1) dimensions for which 
it is possible to obtain the fracture profile of a fibrous 
material submitted to a constant force F, for example, 
by a hanging weight on it. In this work the influence of 
the temperature t on the fracture process was considered. 
We show that, at the boundary between the catastrophic 
and the slowly shredding regime, the roughness exponent 
£ does not depend on the temperature t. 



II. MODEL 

Our model consists of a bundle oi Nq — L x L parallel 
fibers, all with the same elastic constant, k, distributed 
on a triangular lattice. In order to simulate the height of 
the sample, the fibers are divided in r\ segments with the 
same length. The fiber bundle is fixed at both extremes 
by two parallel plates, one of these is fixed and in the 
other a constant force F is applied. This force is equally 
and completely distributed in the fiber bundle, submit- 
ting all fibers to the same linear deformation z = F/Nk, 
where TV is the number of unbroken fibers. This type of 
distribution is called equal load sharing (ELS) P,p|,p3[ . 
Another type is the local load sharing (LLS), where the 
load of a broke fiber is transfer onto its two nearest unbro- 
ken fibers p4] , p5| |. At our model, when the deformation 
z reaches a critical value z e , the failure probability of an 
isolated fiber is equal to one. The failure probability of 
a fiber i is given by 



Pi(6,t) 



1 



• exp 



(S 2 1) 



t 



(1) 



where Ui is the number of unbroken neighboring fibers, 
S = z/z c = F/Nkz c is the strain of the material, 
t = KbT / E c is the normalized temperature, Kb is the 
Boltzmann constant, T is the absolute temperature and 
E c is the critical elastic energy. In this model, besides 
finding the failure probability of a fiber, we have to indi- 
cate in which segment it breaks. Since, each fiber of the 
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bundle can break at different parts a fracture surface is 
produce. Similar procedure was used in Ref to ob- 
tain the fracture profile of a model for fracture on fibrous 
materials in (1+1) dimensions. The segment is randomly 
selected and the probability of the fiber to break in it is 
given by M. 



III. RESULTS 



(j)j{mj 



K + 1) 



(2) 



where uij is a vector which indicates how many times a 
segment j broke and g — J2j( m j + !)• Eq. (||) simulates 
a concentration of tension near to the region where the 
fiber bundle is weaker. 

At the beginning of the simulation, the bundle is sub- 
mitted to an initial strain given by 



S 



F 



N kz c 



(3) 



At each time step we randomly choose a fiber of a set of 



N, 



qN a unbroken fibers. The number q represents a 



percentage of fibers and allow us to work with any system 
size. Then, using Eq. (Q), we evaluate the fiber failure 
probability Pi and compare it with a random number r 
in the interval [0,1). If r < Pi the fiber breaks. We then 
choose a segment j in the fiber and evaluate its prob- 
ability </> to break, using Eq. (||). If the probability <j> 
is higher than a random number / the fiber breaks in 
the chosen segment. If not, we analyze the neighboring 
segments (j+1) and (j-1) and again, evaluate the proba- 
bility </>. If the condition / < <f> does not hold to neither 
of the neighboring segments, we return to the initial seg- 
ment and test the condition / < <f> for a new value of 
/. This process continues until the condition / < <f> is 
true. Once defined the segment where the fiber breaks, 
we begin to test all neighboring unbroken fibers. The 
first segment tested in the neighboring unbroken fibers 
is the one in which the previous fiber broke. The failure 
probability Pi of these neighboring fibers increases due 
to the decreasing of rii and a cascade of breaking fibers 
may begin. This procedure describes the propagation of 
a crack through the fiber bundle in the perpendicular di- 
rection to the applied force. The process of propagation 
stops when the test of the probability does not allow rup- 
ture of any other fiber on the border of the crack or when 
the crack meets another already formed crack. The same 
cascade propagation is attempted by choosing another 
fiber of the set N q . After all the N q fibers have been 
tested, the strain 8 is increased if some fibers have been 
broken. Since the force is fixed, the greater the number of 
broken fibers, the larger is the strain on the intact fibers 
and the higher is their failure probability. Then, another 
set of N q unbroken fibers is chosen and all the rupture 
process is restarted. The simulation terminates when all 
the fibers of the bundle are broken, i.e., when the bundle 
is divided into two parts. 



We performed simulations considering L = 2000 (Nq = 
4 x 10 6 fibers), the elastic constant k = 1, the critical de- 
formation z c = 1 and the number of segments r\ — 1000. 

The failure probability (Eq. ^) can be written as 



r(M) 

(ni + 1) ' 



where the parameter T(t, 8) is defined as 

;8 2 - r 

T{t,6) = Sexp 



t 



(4) 



(5) 



For a triangular lattice (with coordination number 6) 
and r(i, 8) > 6, the rupture of any fiber induces the rup- 
ture of the whole bundle, i.e., the bundle breaks with just 
one crack. Obviously, this crack forms a cluster which 
percolates through the entire system. 

We can define the density of the percolation crack as 



N, 



pc 

TV 



(6) 



where N pc is the number of broken fibers belonging to 
the percolating crack. Thus, when T(t, 8) > 6, we have 
p=l. 

Figure |l| shows the density of the percolation crack 
p versus the initial strain i5o, for two different temper- 
atures. Notice that, for high values of So, p = 1 and 
for low values of Sq the density of the percolating cluster 
p jumps to zero. Thus, we may assume that there is a 
critical value 8q c which depends on the temperature t. 
For the temperatures t = 0.5 and t = 2.0 used in our 
simulation the values obtained for So c are 1.11 and 1.27 
respectively. Above Sq c there is a percolation crack and 
below it any single crack percolates the fiber bundle. The 
critical value, 8q c , represents the transition between two 
failure regimes Q| : catastrophic and slowly shredding. In 
the catastrophic regime there are cracks that percolate 
the fiber bundle and in the slowly shredding any crack 
percolates the bundle. In Ref. |?J we have shown that 
these two regimes are separated by a second order phase 
transition and determined from (Eq. ||) the critical line 
separating these two failure regimes in the plane t x Sq. 

Figure |2| shows the fracture surface obtained for t = 2.0 
and three different initial strains Sq- In Fig. || (a) the 
fracture surface is very rough and this profile is char- 
acteristic of a shredding fracture, in which the speed of 
crack propagation is low due to a slow process of succes- 
sive rupture of fibers in the material. In Fig. || (c) the 
fracture surface presents little roughness and is character- 
istic of catastrophic fracture. Here the crack propagates 
with high speed and the breakage of a single fiber in- 
duces the rupture of the whole the bundle. For Sq = 1.27 
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[Fig. g (b)] the fracture occurs at the boundary between 
the two failure regimes. It can be seen from Fig. || that 
the rupture of the sample begins at different segments, 
since we did not used a deterministic starting notch in 
our simulations. 

In order to evaluate the roughness exponent £ different 
one-dimensional cuts in the fracture surface were consid- 
ered. The roughness W of each cut was found by the 
method of the best linear least-square fitting described 
in [p7| . In this method, the roughness W(e) in the scale 
e is given by 

M 

i=i 

and the local roughness w(e) is defined as 

j=i-e 

a,i(e) and fcj(e) are the linear fitting coefficients to the 
displacement data on the interval [i — e, i + e] centered on 
the fiber i. 

The roughness exponent £ at the critical point <5o c 
for two temperatures t was calculated and a value of 
£ ~ 0.42 ± 0.02 was obtained. Our results indicate that 
at the critical point the value of £ does not depend on 
the temperature t. Figure |^ shows the fits for roughness 
W at t = 0.5 and t = 2.0. We also verified that, as 
the initial strain Sq decreases below 5q c , the rupture pro- 
cess and, consequently, the speed of crack propagation 
become more slow. In this situation the roughness expo- 
nent £ tend to zero. In Fig. ^ we show the plot of the 
time to failure (in Monte Carlo step) Tf as function of 
the initial strain So for t = 2.0. Notice that the time de- 
creases with increase of the initial strain 5q . For 5q > Sq c 
the rupture is catastrophic and the material breaks in the 
first time step. In this regime was not possible to find 
the roughness exponent £. We believe that in this regime 
the fracture surface ceases to be sclf-affinc. 

The roughness exponent £ was calculated along the x 
and y direction. At both directions we verify that the 
roughness exponent has the same value. This behavior is 
expected for a large variety of materials, where the two 
directions have similar scaling properties [^8|,^9| . At the 
transition we have observed a power law behavior on the 
number of cracks versus its size. 

The log-log diagram of the frequency of the cracks H c 
versus their sizes S c observed for t — 2.0 and Sq c — 1-27 
is shown in Fig. S. Note that at the beginning H c seems 
to decay linearly so that we can assume a power law 

H c ~ S~ a , (9) 

where a — 2.032 ±0.007. We observed that at the transi- 
tion this exponent does not depend on the temperature. 



IV. CONCLUSION 



In conclusion, we have studied a model for fracture 
on fibrous materials in (2+1) dimensions which provides 
the fracture surface of the material, in contrast with pre- 
vious models. We calculated the roughness exponent £ 
and showed that it does not depend on the temperature 
t. Our results indicate a value of £ — 0.42 ± 0.02 at the 
boundary between the catastrophic (similar to the dy- 
namic regime) and the slowly shredding regime (similar 
to the quasi-static regime). This value is the same ob- 
tained in other type of simulations in three dimensions 
p0| , pl| that related the exponent £ with the quasi-static 
regime. We also have shown that at the boundary there 
is a power law connecting the frequency H c and the size 
S c of the size. This power law behavior is characteristic 
of systems in the criticality and is a good indication that 
the fracture surface is a fractal. 
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FIG. 1. Density of the percolating cluster p vs the initial 
strain So for two different temperatures: t = 0.5 (circles) and 
t — 2.0 (squares). The data were averaged over 1000 statisti- 
cally independent samples. 
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FIG. 2. Fracture surface for three different initial strains. 
In (a) we have: So = 0.4, in (b) So = 1.27 and in (c) 
So — 1.4. In this particularly simulation we have used a total 
of N = 4 x 10 4 fibers. 
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FIG. 3. Log-Log plot of the roughness W as function of the 
scale e for two different temperatures: t = 0.5 (circles) and 
t — 2.0 (squares). The two straight lines have slope 0.42. 



FIG. 4. Plot of the time to failure (in Monte Carlo step) 
Tf vs the initial strain So for t = 2.0. The data were averaged 
over 1000 statistically independent samples. 





FIG. 5. Log-Log plot of the frequency of the cracks H s vs 
the size of the cracks S c for t = 2.0 and 8o c — 1.27. The 
straight line has a slope 2.032. The data were averaged over 
1000 statistically independent samples. 
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